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ABSmCT 

The  operations  possible  In  a coherent 
optical  computer  are  briefly  reviewed  and  sever- 
al applications  of  this  novel  class  of  processor 
are  Included.  Emphasis  Is  given  to  signal 
processing  operations  and  new  space  variant 
optical  processor  architectures. 


1.  INTRODUCTION 

The  operations  possible  In  a coherent  op- 
tlcal  computer  are  well  known  [1,2]  and  are  thus 
only  briefly  reviewed  In  Section  2.  The  advan- 
* Um  tages  of  optical  computing  Include:  parallel 
li  processing,  high  space  bandwidth,  and  high 
CD  throughput.  The  high  space  or  time  bandwidth 

of  an  optical  system  Is  vividly  demonstrated  In 
1,1  the  folded  spectrim  analyzer  presented  In  Sec- 
^ I tlon  3.  To  demonstrate  the  parallel  processing 
of  these  systems  and  their  high  throughput,  we 
t 1 Include  several  Image  processing  and  signal 
processing  examples  In  Sections  6 and  7.  To 
ai|  realize  the  high  throughput  of  this  processor, 
real  time  and  reusable  spatial  light  modulators 
are  essential.  The  status  of  the  leading  candl- 
date  2-D  transducers  are  briefly  summarized  In 
' ' Section  4.  One  of  the  shortcomings  of  an  opti- 
cal processor  Is  Its  lack  of  flexibility  and  the 
limited  number  of  operations  possible  In  an  opti- 
cal processor.  One  way  of  extending  the  flexi- 
bility of  these  optical  processors  Is  to  combine 
an  optical  and  digital  system  Into  a hybrid  pro- 
cessor. We  briefly  review  current  work  on  such 
systems  In  Section  S.  Another  recent  approach  to 
optical  processing  whereby  the  repertoire  of 
operations  achievable  can  be  Increased  and  many 
practical  problems  overcome.  Is  the  optical  space 
variant  processor.  The  highlights  of  recent  re- 
search In  this  area  are  Included  I.t  Section  8. 

The  multitude  of  topics  and  wealth  of  results  to 
be  discussed  preclude  extensive  description  of 
any  one  application;  these  details  are  available 
In  the  references  noted. 

2,  OPTICAL  COMPUTING 

The  schematic  diagram  of  the  optical  data 
processor  which  we  will  use  to  describe  the  oper- 
ation possible  In  an  optical  computer  Is  shown 


In  Figure  1.  Plane  Pj  Is  the  Input  plane  In  which 
a transparent  version  of  the  Input  data  (signal  or 
Image)  to  be  processed  Is  placed.  Because  we  con- 
sider only  coherent  optical  computers  herein,  the 
spatial  amplitude  transmittance  of  plane  Pj  or  the 
Input  2-D  data  are  assumed  to  be  g(x],yi).  The 
light  distribution  Incident  on  plane  P2  (the  back 
focal  plane  of  lens  Lj)  la  the  2-D  Fourier  trans- 
form of  g or 

G(u,v)  “ //g(Xi,yi)  exp[-j2it  (xju  + yiv)  jdxidyj . (1) 

Capital  letters  (e.g.,  G)  are  used  to  denote  the 
Fourier  transforms  of  the  corresponding  lower  case 
variables  (e.g.,  g) . The  coordinates  (u,v)  of 
plane  P2  are  spatial  frequencies  which  are  rela- 
ted to  the  distance  coordinates  (X2.y2)  of  plane 
P2  by  (X2,y2)  ” (u,v)Xfi,  where  X Is  the  wave- 
length of  the  laser  light  used  and  fj  is  the  focal 
length  of  the  Fourier  transform  lens  Lj. 
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Fig.  1 Schematic  representation  of  a frequency 
plane  correlator. 

In  the  output  plane  P3,  lens  L2  forms  the  > 
Fourier  transform  of  the  light  distribution  G(u,v) 
In  Its  front  focal  plane  F2>  and  In  plane  F3  we 
find  an  image  of  g(xi,yi).  By  placing  simple 
opaque  spots  In  plane  P2,  various  spatial  frequen- 
cies present  In  the  Input  Image  can  be  removed. 
With  an  opaque  spot  In  the  central  or  dc  frequency 
region  of  plane  P2,  the  output  Image  will  be  a 
high-pass  filtered  or  differentiated  version  of 
the  Input  Image.  This  Is  known  as  spatial  filter- 
ing. With  more  complex  spatial  filters  at  plane 
P2  (e.g.,  those  that  spatially  alter  both  the 
phase  and  amplitude  of  the  function).  Images  can 
be  restored  and  enhanced,  blurred  or  defocused 
Images  can  be  restored,  etc. 

With  the  transmittance  of  plane  P2  described 
by  the  Fourier  transform  H(u,v)  of  h(x,y)  and 


r 


with  g(x,y)  placed  at  plane  Fi,  the  light  distri- 
bution Incident  on  plane  P2  Is  G(u,v)  and  the 
light  distribution  transmitted  through  plane  P2  Is 
G«H.  Lens  L2  forms  the  transform  of  this  product 
of  two  transforms  and  at  plane  P3  the  output  pat- 
tern Is  the  convolution  g*h  of  g and  h.  This  Is 
the  basic  operation  used  in  Image  enhancement  and 
restoration  and  Is  equivalent  to  altering  the  Im- 
pulse response  or  point  spread  function  of  the 
system. 

If  we  can  make  the  amplitude  transmittance  of 
plane  P2  proportional  to  the  complex  conjugate  H* 
of  the  Fourier  transform  H of  a reference  function 
h (we  discuss  how  this  Is  possible  below) , then 
with  g(x,y)  at  plane  Pj,  the  light  distribution 
transmitted  through  plane  P2  Is  GH*  and  In  plane 
P3  we  obtain  the  correlation  g^^h  of  the  Input 
and  reference  function.  This  Is  the  basic  opera- 
tion used  In  pattern  recognition  and  In  processing 
coded  waveforms  and  signals  in  radar,  sonar,  etc. 
The  resultant  system  of  Figure  1 Is  then  known  as 
a frequency  plane  correlator  and  the  pattern  H*  at 
plane  P2  as  a matched  spatial  filter. 

To  record  this  H*  pattern  at  plane  P2,  we  use 
the  principle  of  holography.  At  plane  Pi,  we 
place  a transparency  with  amplitude  transmittance 
h(*.y),  l.e.,  the  reference  function.  Its  Fourier 
transform  H(u,v)  produced  by  lens  L]  Is  Incident 
on  plane  P2.  At  plane  P2,  we  record  the  Inter- 
ference of  H(u,v)  and  an  off-axis  plane  wave  refer- 
ence beam  A exp(- j 2ii  0x2 ) (where  A Is  the  uniform 
amplitude  of  the  plane  wave  and  a (sin  6)/!  where 
8 Is  the  off-axis  angle  between  the  reference  bean 
and  reference  function.  The  recording  material  at 
plane  P2  Is  Intensity  sensitive  and  the  subsequent 
transmittance  of  plane  P2  la 

t(x2.y2)  • |H  + A exp(J2irox2)|^ 

" |h|^  + 1a1^  + HA  exp(j2itax2) 

+ H*A  exp  (- j 0X2  ) . (2) 

The  last  term  In  Eq.  (2)  Is  proportional  to  the 
desired  H term  and  can  be  separated  from  the 
other  terms  in  Eq.  (2)  by  proper  choice  of  the 
angle  6 (which  Is  encoded  In  o) . 
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Pig.  2 Schematic  representation  of  a reflex 
readout  Joint  transform  correlator. 


An  alternate  optical  correlator  topology  Is 
shown  In  Figure  2.  In  this  schematic,  the  two 
functions  g and  h to  be  correlated  are  placed 
slde-by-slde  In  plane  P2  separated  by  a center-to- 
center  distance  2b.  The  transmittance  of  Pj  In 
1-D  Is  described  by 

If  (Vl)  “ g(yi-b)  + (yi+b).  (3) 

At  plane  P2  we  record: 

t(v)  - (G(v)  exp(-j2iivb)  + H(v)  exp(j2wvb) 

- |C|^  + |h|^  + GH*  exp(-j4irvb) 

+ G H exp(+j4irvb) . (4) 

Lens  L2  forms  the  Fourier  transform  of  this 
pattern  at  plane  P3,  where  we  find: 

U(x3,y3)  - g*g  + h*h  at  (X3,y3)  - (0,0) 
h at  (X3,y3)  - (0,2b) 
.h©  g at  (X3,y3)  - (0,-2b)  (5) 

where  equal  focal  lengths  fj  - f2  have  been 
asstsned  for  Lx  and  L2. 

Many  versions  of  these  two  systems  exist.  In 
Figure  2,  we  showed  the  recording  of  the  pattern 
at  plane  P2  with  one  wavelength  of  light  and  the 
reading  of  this  pattern  with  a different  wave- 
length of  light  In  transmission  from  beam-splitter 
BS.  The  plane  P2  material  Is  often  read  out  In 
reflection  In  which  case  the  beam  splitter  Is 
placed  to  the  right  of  plane  P2.  In  another 
version  of  either  system,  we  can  record  N signals 
gQ  on  N lines  at  plane  Px  and  replace  lens  Lx 
with  a cylindrical/spherical  lens  pair.  The  pat- 
tern at  plane  P2  then  consists  of  N Fourier  trans- 
forms Gg  of  the  N Input  signals  g^  on  N lines. 

With  a cylindrical/spherical  lens  pei{  for  lens  L2 
and  the  same  N reference  1-D  signal  H on  N lines 
at  plane  F2 , the  output  pattern  at  plane  P3  con- 
sists of  N correlations  g„Wh  on  N separate 
lines.  This  multichannel  1-D  correlator  can  be 
modified  for  use  In  radar  processing  by  use  of  a 
spherical  lens  for  L2.  The  correlation  of  the  N 
input  signals  g^  with  the  reference  signal  h (now 
a coded  reference  radar  waveform)  still  occurs  and 
we  realize  pulse  compression.  The  vertical  trans- 
form down  the  N lines  In  F2  produces  Doppler 
processing  (In  a pulse  burst  radar)  or  provides 
target  azimuth  data  (In  a phased  array  radar) . 
Examples  of  these  and  other  versions  of  these 
basic  correlators  are  Included  In  Sections  6,  7, 
and  8. 


3.  FOLDED  SPECTRUM  OPTICAL  SIGNAL  PROCESSING 
[3,41 

As  one  example  of  optical  signal  processing 
that  vividly  demonstrates  the  large  time  or  space 
bandwidth  product  of  an  optical  processor,  we  con- 
sider a long  1-D  signal  8(t)  of  T sec  duration 
raster  recorded  on  n lines  In  plane  Fx  of  Figure  1 
with  p points  per  line  at  a line  scan  rate 
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f » 1/T 

■•H  H (where  Tjj  is  the  horizontal  line  scan 

time).  The  2-D  Fourier  transform  of  this  pattern 
is  recorded  at  plane  P2  (lens  L2  and  plane  P3  of 
Figure  1 are  not  used  in  this  application).  This 
pattern  is  a folded  spectrum  consisting  of  p hori- 
zontal coarse  frequency  loci  with  n fine  frequency 
points  per  line.  Frequencies  from  0 to  fjj  present 
in  the  input  signal  appear  on  the  first  line  in 
plane  P2  as  spots  of  light,  frequencies  from  fjj  to 
2fjj  on  line  2,  etc.  up  to  a maximum  frequency  pfjj 
on  line  p. 


Fig.  3 Folded  spectrum  output  of  the  RF 

spectrum  with  12  HHz  bandwidth  and 
40  Hz  frequency  resolution  (Courtesy 
Ampex  Corp . ) . 


In  Figure  3,  a folded  spectrum  output  pattern 
of  a 25  msec  record  of  the  RF  spectrum  Is  shown. 
For  this  case,  750  raster  lines  were  recorded  with 
400  cycles/llne.  The  folded  spectrum  pattern  has 
a 12  MHz  bandwidth,  with  400  coarse  Horizontal 
frequency  loci  and  750  fine  frequencies  per  line 
for  a 40  Hz  fine  frequency  resolution  in  a 12  MHz 
signal  bandwidth.  Many  specific  stations  and  fre- 
quencies are  noted  on  the  figure. 


4.  REAL-TIME  REUSABLE  SPATIAL  LIGHT  MODULATORS 

11x61 

Producing  transparencies  of  the  input  g(x,y) 
data  and  the  required  filter  functions  such  as 
H*(u,v)  at  rates  comparable  with  the  throughput 
of  an  optical  processor  requires  real-time  and  re- 
usable spatial  light  modulators  (SLMs)  in  place  of 
photographic  film  at  planes  Pj  and  P2.  The 


desired  data  can  be  recorded  on  these  devices  by 
modulating  the  beam  current  of  an  off-axis  elec- 
tron gun  with  the  Input  signal  in  sync  with  the 
scan  of  the  electron  beam.  These  electron  beam 
addressed  devices  are  useful  for  recording  the 
input  data  at  plane  P^.  The  major  candidate 
electron  beam  addressed  SLMs  and  their  specifica- 
tions are  listed  in  the  top  of  Table  1.  TV  frame 
rates  (33  msec)  are  the  generally  accepted  norm 
for  real  time  although  many  of  these  devices  can 
be  operated  faster. 

The  matched  filter  or  reference  filter  func- 
tion recorded  at  plane  P2  is  usually  easiest 
formed  using  an  optically  addressed  SLM.  These 
devices  generally  consist  of  a photo-conductor  and 
electro-optic  crystal  sandwiched  between  two 
transparent  electrodes  with  a dielectric  mirror 
layer  inside  the  structure.  They  are  generally 
read-out  in  reflection  as  noted  in  the  discussion 
of  Figure  2.  The  leading  candidate  optically 
addressed  SLMs  for  use  at  plane  P2  and  their 
specifications  are  listed  in  the  bottom  of  Table  1. 


5.  HYDRID  OPTICAL/DIGITAL  PROCESSORS  [7,8] 

To  realize  the  high  speed  and  parallel  pro- 
cessing advantages  of  an  optical  processor  as  well 
as  the  flexibility,  programability , and  control 
features  of  a digital  processor,  a hybrid  optical/ 
digital  system  has  been  considered.  The  general 
schematic  of  one  such  system  is  shown  in  Figure  4. 
The  top  portion  of  the  diagram  is  an  optical 
correlator  with  real-time  input  and  filter  plane 
Sl>ls.  The  Fourier  transform  and  correlation 
planes  of  this  system  are  interfaced  to  the  digi- 
tal section  of  the  hybrid  processor  through 
separate  detector  arrays.  The  Interface  digitizes 
the  contents  of  these  planes  into  an  nxm  cellular 
array.  One  function  of  the  rest  of  the  digital 
system  is  to  extract  the  location  and  amplitude 
of  the  "important"  peaks  of  light  in  these  planes 
and  of  course  to  determine  which  are  the  important 
peaks  of  light. 

These  operations  are  achieved  by  thresholding, 
peak  detection,  area  calculations,  interframe 
operations  (addition,  subtraction,  etc.).  These 
are  achieved  by  software  or  hardware  and  a reper- 
toire of  standard  image  processing  operations  in 
Basic.  A PDP-11/15  computer  with  64K  of  core  and 
an  RT-11  operating  system  are  the  master  system. 

An  array  of  32  microprocessors  and  a 4K  250  nsec 
buffer  memory  are  used  to  limit  data  flow  to  the 
PDP-11.  The  digital  preprocessor  consists  of  a 
floppy  disc  system,  a magnetic  tape  unit,  and  a 
2 Mbit  (512  X 512  x 8 bit)  store  whose  total 
contents  can  be  randomly  dumped  in  30  msec.  As 
the  control  lines  in  Figure  4 Indicate,  the  digi- 
tal system  is  used  to  control  the  data  flow,  the 
contents  of  the  SLMs  at  planes  Pj  and  P2,  and  the 
operation  (wrlte/read/or  erase  mode  selection)  of 
the  real-time  SLMs.  The  digital  system  also  con- 
trols the  output  data  displays  (TV,  x-y  CRT, 
isometric  display,  histogram  data,  etc.). 

Examples  of  several  of  these  hybrid  system  opera- 
tions are  included  in  Sections  6 and  7. 
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Table  1.  I^ndldate  real-time  and  reusable  spatial  light  modulators  and  their  specifications. 

The  first  three  are  electron-beam  addressed,  the  last  three  are  optically  addressed. 


E-Beam  DKDP 
TOPR 

Dielectric 

Photo  DKDP 
Prom 

Liquid  Crystal 


RESOLUTION 
1000  X 1000 

1800  X 1800 

500  X 500 

760  X 500 
250  X 250 
650  X 650 


STORAGE 

1 hr. 
days 

30  msec 

2 hrs. 
10  min 
IS  msec 


CYCLE  TIME 

33  msec 
2 sec 
33  msec 

1 msec 
1 msec 
33  msec 


BANDWIDTH 

20  MHz 
70  MHz 
10  MHz 
SENSITIVITY 
10  ergs/cm^ 
5000  ergs/cm^ 
5 ergs/cm^ 


SHUTTER,  P, 


OISITAL 

prerrocessor 


ARRAY  a BUFFER 


radiation  pattern  are  Fourier  transform  pairs. 

The  dual  relationship  states  that  the  received 
signal  pattern  and  the  location  of  the  radiating 
object  are  also  Fourier  transform  pairs.  For 
simplicity,  consider  a linear  phased  array  of  N 
elements.  If  the  time  history  of  the  received 
signals  from  the  N elements  of  the  array  are 
recorded  on  N lines  on  the  electron-beam  KDKF  at 
plane  P^,  the  2-D  Fourier  transform  of  this  Input 
pattern  (appearing  at  plane  P2)  consists  of  a 
central  dc  term  and  two  off-axla  peaks  of  light. 

A typical  output  Fourier  transform  pattern  Is 
shorn  In  Figure  5a.  The  horizontal  coordinates  of 
the  t«ro  off-axis  peaks  of  light  are  proportional 
to  the  target’s  azimuth  angle  and  thus  represent 
the  desired  output  plane  data. 

After  digital  processing  to  remove  the 
central  dc  term  and  determine  the  coordinates  of 
the  off-axis  terms  In  Figure  Sa,  the  simpler  dis- 
play of  Figure  Sb  results.  Extensive  phased  array 
data  was  processed  and  In  all  cases  excellent 
agreement  with  the  predicted  values  obtained  (see 
Figure  Sc).  The  optical  processor  easily  accom- 
modates multiple  targets.  In  this  case,  multiple 
output  peaks  of  light  result;  an  example  of  this 
for  pulsed  Doppler  radar  follows.  Planar  phased 
array  radar  data  has  been  processed  by  similar 
methods;  In  this  case,  the  coordinates  of  the  out- 
put off-axis  peak  of  light  are  proportional  to  the 
target's  elevation  and  azimuth  angles. 


6.2  Pulsed  Radar  Optical  Signal  Processing  fill 

To  optically  process  pulsed  Doppler  or  FM 
stepped  radar  data,  N received  signals  are  written 
on  N lines  on  the  electron  beam  DKDP  device  at 
plane  Pj  (with  a horizontal  scan  rate  equal  to  the 
radar's  prf)  and  the  2-D  Fourier  transform  of  this 
Input  pattern  viewed  at  plane  P2.  The  output 
pattern  Is  similar  to  that  of  Figure  Sa  except  that 
the  horizontal  coordinates  of  the  off-axis  peaks  of 
light  are  now  proportional  to  the  target’s  fine 
I velocity  (Doppler  radar)  or  fine  range  (FM  stepped 

radar) . In  Figure  6a,  the  output  plane  pattern 
for  real  pulsed  Doppler  radar  data  Is  shown.  The 
two  pairs  of  off-axis  peaks  of  light  correspond  to 
two  aircraft  crossing  In  the  same  range  bln  at  two 
different  Doppler  velocities.  This  pattern,  after 
digital  processing.  Is  shown  in  Figure  6b  from 
which  the  velocities  of  the  two  targets  have  easily 
been  found. 


Fig.  6 Real-time  pulsed  Doppler  radar  data 

optical  processing  using  electron  beam 
DKDP.  The  horizontal  coordinates  of  the 
off-axis  spots  are  proportional  to  the 
target's  Doppler  velocity,  (a)  Optical 
output  pattern;  (b)  Digitized  version  of 
(a).  [11] 


Pj.  The  2-D  transform  of  the  light  distribution 
transmitted  through  plane  P2  now  has  the  compres- 
sion and  high  slgnal-to-noise  of  matched  filtering 
and  coded  waveforms.  The  coordinates  of  the  out- 
put peak(s)  of  light  are  now  proportional  to  the 
target's  azimuth  and  range  (for  a phased  array)  or 
the  target's  Doppler  and  range  (for  a pulse  burst 
radar)  and  thus  constitutes  the  ambiguity  function 
of  the  radar.  Two  examples  are  shown  in  Figure  7 
for  a target  at  azimuth  angles  of  5°  and  30°  and 
fine  ranges  of  45%  and  55%  of  the  range  bln  Rg. 


6. 3 Range/Doppler  or  Azimuth/Range  Optical 
Signal  Processing  [121 

To  extend  this  optical  signal  processing  tech- 
nique for  use  with  coded  radar  waveforms,  the  full 
system  of  Figure  4 Is  used  with  a cylindrical/ 
spherical  lens  pair  for  Lj . N Identical  reference 
signals  of  ths  coded  radar  waveform  are  recorded 
on  N lines  at  Pi  and  a multi-channel  1-D  matched 
spatial  filter  of  the  coded  reference  waveform 
recorded  at  plane  P2.  N pulse-burst  radar  returns 
or  the  received  signals  from  N elements  of  a phased 
array  are  then  recorded  on  N lines  at  plane 
t 


6.4  Optical  Correlation  of  Long  Coded  Waveforms 

[13]  A 

The  radar  code  Is  often  quite  long  (e.g.  10 
bits)  and  the  linear  resolution  of  an  SIM  cannot 
accommodate  it.  Another  problem  associated  with 
such  a code  is  the  determination  of  where  within 
the  transmitted  code  a received  code  starts.  For 
a long  NM  bit  code,  we  devised  an  optical  proces- 
sor capable  of  solving  both  of  the  problems.  The 
received  signal  Is  raster  recorded  on  M lines 
with  N points  per  line.  A unique  matched  filter 
of  only  (aj-l)x(2N-l)  bits  was  constructed.  This 
Is  far  less  than  the  (NMl^  bit  size  normally 


Multi-channel  1-D  optical  processing 
for  range/azimuth  data  (a)  Target  at 
5“  and  45%  of  the  range  bln;  (b)  Target 
at  30°  and  55%  of  the  range  bln.  [12] 


required.)  The  system  Is  most  useful  for  pseudo- 
random codes.  The  output  correlation  plane  pattern 
consists  of  a spot  of  light  of  full  amplitude  with 
Its  location  proportional  to  the  starting  bit  lo- 
cation of  the  received  code  within  the  full  trans- 
mitted code. 


7.  OPTICAL  IMAGE  PROCESSING  [7,141 


In  this  brief  section,  we  have  selected  sev- 
eral image  processing  operations  that  can  be  per- 
formed directly  on  a stored  image  using  photo  DKDP 
and  two  optical  pattern  recognition  or  correlation 
examples.  These  latter  examples  were  again  per- 
formed on  the  hybrid  processor  of  Figure  4 in  real 
time  using  the  electron-beam  addressed  DKDP  SLM  in 
the  input  plane  Pj. 


7.1  Image  Processing  [141 


Three  image  processing  operations  performed 
in  real  time  on  an  Image  stored  on  the  photo-DKDP 
SLM  without  the  need  for  external  spatial  filters 
etc.  are  described.  In  Figure  8a,  we  show  an 
example  of  image  addition.  The  two  images  were 
focussed  onto  the  device  using  two  separate  expo- 
sures with  the  same  polarity  of  control  voltage 
across  the  DKDP.  With  this  multiple  Image  pattern 
stored  on  the  photo-DKDP,  the  SLM  was  exposed  to 


the  second  image  with  the  polarity  of  the  voltage 
across  the  photo-DKDP  reversed.  This  results  in 
the  subtraction  of  two  Images  as  shown  In  Figure 
8b.  If  the  control  voltage  across  the  SLM  is 
varied  while  an  image  (Figure  8c)  is  present,  any 
image  contours  present  in  the  image  can  be  dis- 
played shown  in  Figure  8d. 


Fig.  8 Real-time  image  processing  using  photo- 
DKDP.  (a)  Image  addition;  (b)  Image 
subtraction;  (c)  Input  image;  (d)  Con- 
tour generation  on  the  Image  In  (c)(14] 
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7.2  Optical  Pattern  Recognition  [71 

With  an  aerial  Image  (Figure  9a)  recorded  on 
the  electron-beam  DKOP  at  plane  Pj  of  Figure  4 and 
a matched  spatial  filter  of  the  desired  reference 
function  (Figure  9b)  recorded  at  plane  P2,  the 
output  pattern  at  plane  F3  Is  the  correlation  of 
the  Input  and  reference  functions.  The  amplitude 
of  the  output  correlation  peak  Is  proportional  to 
the  agreement  between  the  Input  and  reference 
functions.  The  location  of  the  output  correlation 
peak  is  proportional  to  the  location  of  the 
reference  function  In  the  Input  scene  as  shown  In 
Figure  9c.  The  binary /digitized  version  of  this 
pattern  (Figure  9d)  Is  easily  determined  by  the 
digital  system.  The  location  of  the  correlation 
peak  can  be  used  as  feedback  control  In  a missile 
guidance  application  In  which  the  Input  image 
originates  from  an  on-board  sensor  and  the  refer- 
ence Image  is  one  from  a large  stored  data  base. 


Fig.  9 Real-time  optical  pattern  recognition 
of  aerial  Imagery  using  electron-beam 
DROP,  (a)  Input  Image;  (b)  Reference 
Image;  (c)  Output  correlation  pattern; 
(d)  Digitized  version  of  (c)  [7] 


As  noted  earlier,  the  optical  correlator  of 
Figure  4 can  be  used  with  multiple  occurences  of  a 
reference  function.  As  an  example  of  this,  the 
Input  pattern  of  Figure  10a  was  recorded  on  the 
electron-beam  DKDP  at  plane  F;  of  Figure  4 and  a 
matched  spatial  filter  of  the  word  "Cod"  formed  at 
plane  P2 . The  output  correlation  pattern  at  plane 
P3  Is  shown  In  Figure  10b.  After  digital  process- 
ing, the  desired  pattern  of  Figure  10c  results  In 
which  the  cross-correlation  typical  of  text  recog- 
nition have  been  removed  and  where  the  locations 
of  each  of  the  seven  correlation  peaks  are  propor- 
tional to  the  locations  of  the  seven  occurrences 
of  the  word  "God"  In  the  Input  Image. 


Fig.  10  Real-time  optical  pattern  recognition  of 
words  using  electron-beam  DKDP.  (a) 
Input  pattern;  (b)  Optical  correlation 
pattern  of  the  word  "God"  and  the  image 
In  (a);  (c)  Digitized  version  of  (b). 

[7] 


8.  SPACE  VARIANT  OPTICAL  PROCESSING  [15,161 

To  increase  the  flexibility  of  an  optical 
processor  and  avoid  several  sources  of  correlation 
degradation  such  as  scale  and  rotational  changes 
between  the  input  and  reference  functions,  a 
space  variant  optical  correlator  has  been  sug- 
gested. A coordinate  transformation  gjfx)  ■+ 
gl(e*)  Is  applied  to  the  input  function  (1-D  func- 
tions are  used  for  simplicity  only)  (here  a simple 
log-scaling  of  the  coordinates  is  used),  and  this 
new  gi(e*)  function  is  recorded  at  Pj  of  Figure  1. 
The  Fourier  transform  of  this  function  formed  at 
P2  Is  the  Mellln  transform  Mi(u)  of  gi(x).  With 
M;  recorded  at  plane  P2 , and  a scaled  version 
g2(ax)  of  gi  recorded  at  P;  as  g2(e®’'),  the  light 
distribution  transmitted  through  P2  is  Mi*M2  where 
M2  ” aJ'*Mi  “ M;  exp(Julna),  Lens  L2  forms  the 
Fourier  transform  of  Mi*M2  and  at  plane  P3  we  find 

UsCx)  “ 81  © 82  “ ®1  © ®1  * (6) 

From  Eq.  (6),  we  see  the  power  of  this 
syst^.  The  magnitude  of  the  output  correlation 
gl  ® g2  of  two  functions  that  differ  in  scale  is 
Identical  to  the  auto-correlation  gi®  g]  of  gj 
and  thus  there  is  no  loss  in  signal-to-noise  ratio 
and  we  have  produced  a scale  Invariant  correlator. 
The  delta  function  In  Eq.  (6)  indicates  that  the 
location  of  the  correlation  peak  is  proportional 
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1 


1 
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to  the  scale  change  between  the  Input  and  refer- 
ence functions. 


t 
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Fig.  11  Scale  and  rotation  Invariant  correlation 
(a)  Input;  (b)  Reference;  (c)  Output 
correlation;  (d)  Cross-section  of  auto- 
correlation; (e)  Cross-section  of  cross- 
correlation  tl5] 


between  them,  we  can  Mellin  correlate  these  sig- 
nals with  no  loss  in  slgnal-to-nolse  and  from  the 
location  of  the  correlation  peak  determine  the 
target's  Doppler.  An  example  of  this  using  Che 
correlator  of  Figure  2 with  cylindrical/spherical 
lens  pairs  for  Lj  and  1,2  shown  in  Figure  12. 

The  same  nine  signals  were  recorded  on  one  half  of 
the  input  plane  and  nine  different  Doppler  ver- 
sions of  the  signal  on  the  other  half  of  the  in- 
put plane  (with  a 100%  Doppler  difference  across 
the  signals) . The  output  correlation  plane 
pattern  of  Figure  12  shows  nine  correlation  peaks 
of  equal  intensity,  with  the  horizontal  position 
of  each  spot  proportional  to  the  Doppler  differ- 
ence between  the  Input  signal  on  the  respective 
line  and  the  reference  waveform. 


i 
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i 

I 
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A log-coordinate  transform  results  in  a scale 
Invariant  correlator  as  described  above.  To 
realize  a rotational  invariant  correlator,  a polar 
coordinate  transform  f(x,y)  -*•  f(r,9)  is  performed. 
With  p “ e*^  and  the  coordinate  transformation 
f{*«y)  f(p,6),  a scale  and  rotation  Invariant 

correlator  results  where  the  slgnal-to-nolse  of  the 
output  correlation  Is  Invariant  to  scale  or  rota- 
tion changes  between  the  Input  and  reference  func- 
tions and  where  the  coordinates  of  the  output  cor- 
relation peak  are  proportional  respectively  to  the 
scale  change  and  rotation  differences  between  the 
Input  and  reference.  The  correlation  of  the  two 
inputs  (Figures  11a  and  b)  that  differ  In  scale  by 
200%  and  are  rotated  by  e,,  - 180”  Is  shown  In 
Figure  11c  and  the  cross-sectional  scans  of  the 
auto-correlatlon  of  one  object  and  the  cross- 
correlation of  both  objects  in  Figures  lid  and  He. 
As  shown  there  Is  no  slgnal-to-nolse  ratio  loss  and 
the  displacement  of  the  location  of  the  correlation 
peak  Is  proportional  to  the  scale  change  between 
the  objects. 

The  Mellin  transform  has  Immense  practical  use 
In  Doppler  signal  processing.  Since  two  signals 
that  differ  in  Doppler  have  only  a scale  change 


Fig.  12  Optical  Doppler  signal  processing  using 
Mellin  transforms.  [16] 


9.  SUMMARY  AMD  CONCLUSIOM 

From  the  selected  examples  above,  the  high 
space  bandwidth,  parallel  processing,  and  signal 
and  image  processing  applications  of  an  optical 
processor  are  apparent.  The  hybrid  and  space 
variant  processors  described  offer  great  promise 
for  such  systems  and  the  development  of  novel 
data  processing  schemes. 
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